January 9, 2003 



The Gribov problem and its solution from 
a toy model point of view 

Kurt Langfeld 



Institut fur Theoretische Physik, Universitat Tubingen 
D-72076 Tubingen, Germany. 



Abstract 

The standard Faddeev Popov gauge fixing procedure is put on solid grounds by 
taking into account a topological factor which corrects for the number of Gribov 
copies within the first Gribov horizon. Zwanziger's stochastic approach to gauge 
fixed Yang-Mills theory is briefly reviewed. A simple toy model is presented which 
illustrates both methods. Within the toy model, I show that a stochastic drift force 
can be constructed with which the gauged configurations are attracted by the fun- 
damental modular region. The toy model shows that an action which gives rise to 
the drift force can be found. This makes a "heat bath" simulation possible, which is 
seen to be superior to the Langevin approach at the numerical level. 
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1 Introduction 



Although lattice gauge simulations can directly address gauge invariant, physical observ- 
ables of the theory of strong interactions, it might turn out fruitful for detecting basic 
mechanisms to remove redundant degrees of freedom by means of gauge fixing. Moreover, 
in the case of the perturbative approach to QCD Greenf unctions as well as in the case 
employing Dyson- Schwinger techniques 5112], gauge fixing is inevitable. 

In the Faddeev Popov approach, gauge fixing is imposed by demanding that the gauge 
fields obey the gauge fixing condition. The degeneracy factor which determines the weight 
with which a representative of the gauge orbit contributes to the partition function is 
assumed to be given by the Faddeev Popov determinant. It was firstly pointed out that 
e.g. the well known Landau gauge condition is not sufficient for an unambiguous choice 
of the gauge field 3J. As a further refinement of the gauge condition, one might demand 
that the lowest eigenvalue of the Faddeev Popov matrix evaluated with gauged fields, is 
positive [3 E]. These gauge fixed configurations are said to lie within the first Gribov 
horizon. It still turns out that more than one representative of the gauge orbit lies within 
the first Gribov regime in the generic (non-perturbative) case 0. Finally, the domain of 
gauged configurations which have been unambiguously selected from the first Gribov regime 
is known as the fundamental modular region. In practice, any of the refinements of the naive 
Landau gauge condition can be used for a calculation of the gauge invariant observables, 
once the weight factor of each configuration is known. The Gribov problem emerges when 
this weight factor is calculated by the standard Faddeev Popov method: it was pointed 
out in [H1IIIIH! that the inverse Faddeev Popov determinant identically vanishes leaving us 
with an ill-defined partition function. There are two possibilities which are actually used 
in lattice gauge theories to avoid this problem: (i) gauge fixed configurations are randomly 
selected from the first Gribov regime, and their weight factor is exactly taken into account 
(see e.g. [H]-[I2]); this method does not suffer from the Faddeev-Popov problem, but the 
selection of the configurations from the first Gribov regime is still arbitrary. This could 
lead to ambiguous results for gauge variant quantities, such as the gluon propagator; (ii) 
a Laplacian version [T3] of the gauge condition is chosen which allows a unique definition 
of the gauged configuration in practical simulations (see e.g. [13]). 

In order to avoid the Gribov problem in an ab initio continuum formulation of Yang-Mills 
theory, it was proposed by Zwanziger |U [T3] to select stochastically, but in a well defined 
manner the gauge field configurations from the configuration space. Thereby, preference is 
given to configurations within or close to the first Gribov regime. Using the framework of 
the stochastic quantization ^H], the bias towards the first Gribov regime is provided by a 
drift force "tangent" to the gauge orbit. 

In the present paper, I briefly review the gauge fixing techniques of continuum Yang-Mills 
theory. A simple toy model is used to illustrate the Gribov problem of the standard Faddeev 
Popov quantization, and its resolution in the stochastic approach. I find in the case of the 
toy model that a modification of the stochastic approach generates configurations which 
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are attracted by the fundamental modular region. I finally point out that a "potential" 
representation of the drift force seems possible. This paves the path to a very efficient 
numerical simulation using the heat bath techniques. 



2 Brief review of gauge fixing 

2.1 The standard gauge fixing procedure 

The task boils down to select a single gauge field configuration of {A^(x)}, where 
is the set of gauge fields A^(x) which have been generated from a representative A^x) 
of the gauge orbit by applying all possible gauge transformations Q(x). Usually, a gauge 
fixing condition, e.g. such as Landau gauge 

fi(aO : d.A^x) = , (1) 

is chosen to accomplish this task. Unfortunately, this procedure is plagued by the so- 
called Gribov problem [3]: background fields A^x) which give rise to "large" field strength 
generically admit several solutions to the gauge fixing condition (for an illustration see [TT?]). 
This implies that the topological field theory, which generates to the inverse Faddeev- Popov 
determinant, vanishes identically [HI 13 El 

i _ / T->n x ( a 



A-/ P = J Vn5[d^{x)) = 0. (2) 

In the present paper, a lattice regularization as well as a finite space-time volume is im- 
plicitly understood in order to render functional integrals well defined. In order to select 
a gauge configuration A^(x) from the possible solutions of (JTJ), one might demand 



f[A n ] := <V» = , M := , A > , (3) 



where Ao is lowest eigenvalue of the operator M.. The configurations A^(x) emerging from 
the latter equation are said to be restricted to the first Gribov regime. It turns out, that 
even the refinement of gauge fixing does not yield a unique representative A^(x) of 
the gauge orbit. In order to count the number n[A] of configurations lying within the first 
Gribov horizon, one evaluates the functional integral, i.e., 



n[A] = JvQ 5{d,A^x)) Det(^l) 9 



Ao 



where the modulus of the functional determinant has been replaced by the determinant in 
view of the constraint to positive eigenvalues only. Since the gauge condition (jSJ) can be 
installed using the variational condition 



d A x A 2 ^(x) — > minimum 
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it is clear that at least one representative can be found on the gauge orbit implying that 
n[A] > 1. Hence, the partition function of gauged fixed Yang-Mills theory is well defined 
and is given by 

Z = JvAn-^S^A^ Det(^l) #[a ] exp{-S YM } . (4) 

The determinant in the latter equation can be represented by ghost fields as usual. At 
the current stage of investigations present in the literature the topological factor n[A] is 
set to unity. A perturbative treatment of the partition function is consistent with the 
approximation n[A] = 1, since only "small" field strengths are involved by construction. 
The influence of the factor n[A] on a genuine non-perturbative approach is not known. To 
my knowledge, an analytic expression of the integer n[A] is not available. 

In order to circumvent these difficulties, Baulieu and Schaden proposed [H|-[E] to perform 
subsequent steps of gauge fixing where each of them avoids the topological obstruction ©. 



2.2 The stochastic approach to gauge fixed YM-theory 



It was pointed by Zwanziger |HEI] that the Gribov problem can be avoided when stochastic 
quantization ^B] is used in a modified form: a Langevin simulation generates a series 
of configurations {A^(x)} L where each configuration possesses a certain bias concerning 
its position on the gauge orbit. At the same time, the correct result for gauge invariant 
observables is recovered when the observable is calculated from the Langevin "time" series. 
In order to be specific, let me firstly point out that Langevin time series is generated by 



9 , / x 



5S- 



YM 



SA^x) 



+ JC^(x) + [noise] 



(5) 



where Sym is the Yang-Mills action, and fC^x) is a drift force, which will be specified 
later. Expectations values and Greenfunctions are evaluated using ^H] 









A^x; t) 


Jo 



VA,Q[A]P[A] cc hm - 
The probabilistic weight P[A] satisfies a Fokker-Planck equation 



H FP P[A] 



<Px 



S 



SA^x) 



5P[A] 



SA^x) 



+ - 



5S 



YM 



SA^x) 



P[A] 



(6) 



0. (7) 



For a vanishing Drift force, i.e., /C M (x) = 0, one recovers the desired result 



P[A] oc expj-SViu} 



The key observation ^3] made by Zwanziger is that if one chooses a drift force "tangent" 
to the gauge orbit, the gauge fields are generated with preference to a certain region of the 
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gauge orbit, and, hence, are interpreted as gauge fields of a certain gauge. In the context 
of Landau gauge, one chooses 

/C M = D^v(x) , 

where is the gauge covariant derivative and v (x) is an arbitrary auxiliary field On 
the other hand, gauge invariant observables are independent of /C M (x), and, therefore, the 
right hand side of © reproduces the correct Yang-Mills result for gauge invariant quanti- 
ties. The goal of the work ^1 is that a theoretical, Dyson- Schwinger type framework is 
developed for the stochastic approach. A numerical analysis using the stochastic approach 
can be found in jT3 HE] . 

3 The toy model 

3.1 Settings 

Instead of a functional integral over the gauge fields A^x), I will study a simple integral 
of two variables x\ = x, X2 = y. Gauge invariance is replaced by the rotational invariance 
in two dimensions. The "gauge invariant" action is given by 

Sym = x 2 + y 2 = x\ . (8) 

In this very simple toy model, a separation of the degrees of freedom into gauge invariant 
and gauge dependent ones is done by introducing polar coordinates (r, ip) . "Gauge invariant 
observables" are functions of r only. We are interested in expectation values of "gauge 
invariant" observables: 

(f(r)) = ^ J dxdyf(r) exp{-S YM } , (9) 

where r := ^ x 2 + y 2 . If x denotes a representative of a gauge orbit, the orbit {x^} is 
generated by 

/ COS6 — Sm# \ _ r r - . 

x u = ) x , 6 e 0,2tt . 10 

V smu cos v I 

3.2 Standard gauge fixing 

In order to remove the rotational degree of freedom, we introduce a gauge fixing condition 

y(9) = . (11) 
One easily verifies that the procedure suffers from the topological obstruction 0, i.e., 

A^p = J d6 5(y{0)} = J d6 5 (sin6 x + cos6y ^ = (12) 
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independent of the gauge orbit specified by (xo,yo)- To simplify the discussion below, we 
choose Xq > 0, y = without a loss of generality. In order to refine the gauge fixing 
procedure according to (J!J), we calculate 



From 



dy{9) 



dy(0) 
dd 



cosOxq . 



(13) 



de 



COS 9 Xq > 



we conclude that the first Gribov regime consists of the positive half-space x > 0. Conse- 
quently, the topological quantity n[A] becomes 



n[A] 



n\x\ 



it/2 



dO cos 9 5 ( sin xq 



it/2 



1 . 



(14) 



In the case of the example ([lip, there is only one configuration within the first Gribov 
horizon. In this case, the first Gribov regime coincides with the fundamental modular 
region. 



Let us study the more advanced example of a gauge fixing condition: 

1 



f(9) = x(9)y(9) = 0, 



Xq sin# cos 9 



Xq sin(2#) = . 



In this case, four gauge copies exist. The first Gribov regime is determined from 

df(9) 



dO 



2 cos(2#) > 



9 G [0,2vr[ . 



(15) 



(16) 



The space where (j!6|) holds is shaded in figure H Hence, the first Gribov regime consists of 
the complete x-axis. A direct calculation of n[x] shows that there are still two gauge copies 
within the first Gribov horizon (see figure CJ. Nevertheless, we arrived at a stage where 
calculations of "gauge invariant" quantities are feasible: using n[x] =2 and adapting (j3J), 
we find _ 

x 2 — y 2 exp< —x 2 — y 2 \ . . (17) 



Z = 2tt J dx dy - Syxyj 9 

The result that n[x] is independent of Xq is due to an oversimplification by the toy model. 
One can easily construct examples for the generic case, i.e., n[x] is a function of Xq. 



3.3 The stochastic approach 



In the present case of only two degrees of freedom, the gauge invariant part of the Fokker- 
Planck Hamiltonian (|7J) is given by a partial differential equation 



Hir 



d_ 

dx,. 



d_ 

dx,. 



dS 



YM 



dx,. 



(1* 
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Figure 1: Illustrations of the configuration space using the "gauge condition" (fI5jl. 



Infinitesimal rotations are generated by the operator 





G 



dx k 



(19) 



where is the total anti-symmetric tensor in 2-dimensions. It is easy to check that the 
operator G (^TQj) commutes with H inv . Following the procedure suggested by Zwanziger [TK] . 
the drift force "tangent to the gauge orbit" is chosen as Ki = €i k x k v(x), where the scalar 
function v{x) will be specified below. The Fokker-Planck equation (which determines the 
weight factor P v ) with the restoring drift force included is given by 



HfpPv 



d_ 

dxi 



dxi 



+ 



dS 



YM 



dxi 



+ e ik x k v (x) )P V 



. 



(20) 



Case study I: The trivial case. 

Let us study the "gauge invariant" scalar function v(x) = a _1 /r 2 . Using this particular 
choice, the Fokker-Planck equation is given by 



-d 2 - 2d k x k + [x 2 d x -x 1 d 2 



. 



To shed light onto the solution of (|2ip. we introduce polar coordinates: 
x\ — r cos ip , x 2 = r sin , G [0, 2 r w[ . 

One finds: 

1 



-d 2 



i_ L 2 _ _ 4 + _^ L3 

r z a r z 



(21) 



(22) 



(23) 
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Figure 2: The expectation values (r 2 ) and (y 2 ) as function of the Landau parameter a -1 
(left panel). The probability distribution of x n for a -1 = 5 (right panel). Landau type 
drift force, see (|2*7|) . 

where 

L 2 = d% L 3 = . (24) 

As in quantum mechanics, the wave function P v factorizes (for this particular choice of 
v(r)): 

P v = e imLp P v (r), m 2 + iaT^rn = 0. (25) 

The only solution obeying periodic boundary conditions, i.e., P v ((p + 2n) = P v (ip), is m = 0. 
The Fokker-Planck solution is hence given by 

P v = exp|-r 2 | . (26) 

One recovers the trivial result of the "un-fixed" action. 
Case study II. The Landau gauge. 

Let us now investigate a drift force which is not constant along the gauge orbit, i.e., the 
"Landau type gauge" is given by the condition x\x e 2 = where the vector x e k is generated 
from the vector x^ by a rotation by an angle 9 (see (jll))) ). Here we choose the scalar 
function of the drift force according v(x,y) = a^ 1 xy, i.e., 

Ki = e ik x k v(x) , v(x)=a~ 1 xy. (27) 

The direction of this drift force tangent to the gauge orbit is illustrated in figure ^ (right 
panel). 

Rather than seeking the solution of the Fokker Planck equation for the case (|27|). I employ 

(n) 

the corresponding Langevin equation to generate a "time" history x k of the vectors and 
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to calculate expectation values. The "time" history is obtained by the recursion 

x (n+i) = x (n) + dt f-?^L + a^x/j + V, , (28) 

V {n+l) = V {n) + dt (-^ - a-** 2 *) + Vy , (29) 

where r] x / y is a Gaussian noise with width yAdt. An extrapolation dt — > must be 
performed when expectations values are calculated. For the present example, N = 10 7 
pairs (x n ,y n ) are created for dti = 0.01 and dt 2 = 0.005, respectively. The estimator, 
s(dt), of a desired observable is calculated and the corresponding statistical error Ss(dt) is 
estimated. The final error e comprises statistical as well as systematic errors due to the 
extrapolation dt — > 0. I used the following procedure: 

s = -(s{dh) + s(dt 2 )^ , Ss 2 = 5s 2 (dt l ) + 5s 2 (dt 2 ) + {s{dt 1 ) - s(dt 2 )) 2 . 

For a smaller value of dt, one must choose a larger value of iV of configurations in order to 
explore the complete "configuration space". Generically, the Langevin approach is by far 
inferior as e.g. the standard heat bath approach. 

Let us firstly study the "gauge invariant" operator (r 2 )(a~ l ) as function of the Landau 
parameter a. The numerical result is shown in figure El The error bars comprise statistical 
as well as the systematic error from the extrapolation dt — > 0. One indeed finds that (r 2 ) 
is independent of a within numerical accuracy. I also checked that (r 4 ) is independent of 
of a. 

Secondly, I study the "gauge variant" expectation value (y 2 ). The result is also shown in 
the above figure. One finds that (y 2 ) = (r 2 )/2 for a" 1 = as it should be. Increasing 
a -1 decreases fluctuations around y = 0: the configurations (x,y) are pushed to the first 
Gribov regime specified by y = 0. I then calculated the expectation value of (x). It turns 
out that it vanishes for the investigated range of a G [0, 10]. This shows that there is still an 
average over the complete first Gribov regime, which consists of the positive and negative 
half x-axis. The probability distribution of x is also shown in figure 121 for a" 1 = 5. The 
values x^ n > are symmetrically distributed around x\ — x — 0. Hence, the configurations 
are sampled over the complete first Gribov regime (rather than the fundamental modular 
region, i.e., the positive half x-axis). 

Case study III. The fundamental modular region. 

In this section, I will study a drift force which pushes the configurations (x n ,y n ) towards 
the fundamental modular region y = 0, x > 0. The drift force is given by 

Ki = e ik x k v(x) , v (x) = -a' 1 y/r 3 , r = a/x 2 + y 2 . (30) 

The numerical simulation reveals that the "gauge invariant" observable (r 2 ) is independent 
of the gauge fixing parameter a" 1 . In addition, the points (x^ n ',y^) are nicely attracted 
by the positive half x-axis (see figure EJ). 
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Figure 3: The expectation values (r 2 ) and (y 2 ) as function of the Landau parameter a 1 
(left panel). The probability distribution of x n for a -1 = 5 (right panel). Drift force in 



The crucial observation in deriving the function P v (x) is that the drift force ()3())1 can be 
written as the gradient of a gauge fixing action function, i.e., 



Kt — — diSfa , Sfix 
Hence, the Fokker-Planck equation is given by 



a 1 x/r . 



HfpPv 



_d_ 

dxi 



dP v 



+ 



dS 



YM 



<9£ 



fix 



Pv 



dxi \ dxi dxi 
From this representation, one can read of the desired solution, i.e. 

P v (x) = expj- 



-r 2 + a 1 cosip 



(31) 



(32) 



(33) 



Here, the angular part and the radial part factorizes. This proves that gauge invariant 
observables are indeed independent of the gauge fixing parameter. The maximum proba- 
bility of P v (x) is obtained for ip = 0, which is the positive half x-axis, i.e., the fundamental 
modular region. 

The heat-bath simulation. 

Finding the gauge fixing action S^ x which generates the drift force paves the path to 
an efficient simulation using the heat-bath approach. Thereby, configurations (x, y) are 
generated according the probability 



expj-SVM - 5 fix | 



(34) 
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Figure 4: The expectation values (r 2 ) and (y 2 ) as function of the Landau parameter a -1 
from a heat-bath simulation (left panel). The scatter plot of 10 4 heat-bath configurations 
for a -1 = 5 (right panel); Drift force in (|30|) . 

where Sfi x is given by (p?Tj) . It turns out that N h = 10 4 heat-bath configurations are 
sufficient for error bars at the 1% level (compare figure El with figure EJ). A scatter plot of 
the configurations (x n ,y n ) is also shown in figure (JHJ). 

In conclusion, the heat-bath approach is roughly two orders of magnitude more efficient 
than the corresponding Langevin-approach. It is therefore highly desirable to construct 
the gauge fixing function Sq x generating the drift force. 

For these purposes, the drift force Ki must obey certain constraints in order to be generated 
by a gauge fixing action. These constraints are derived from (|31|) by observing that 

- (d k di - did^j = , (35) 

d k K { - diK k = . (36) 

The only non-trivial information is obtained from (}3l)j) by choosing i = 1 and k = 2. Using 
Ki from (|3(jp . we find 

xd x v(x,y) + yd y v(x,y) + 2v(x,y) = 0. (37) 

The last expression is most instructive using polar coordinates, and ()37|) becomes 

r d r In v(r,(p) + 2 = 0. (38) 

Hence, the general solution of v(r,(p) which admit a representation in terms of a gauge 
fixing function is given by 

v(x,y) = v(r,<p) = f(<p) , (39) 
where the function f{<p) is still arbitrary. 
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4 Conclusions 



A simple toy model was designed to illustrate the Gribov problem of the standard Fad- 
deev Popov quantization. Using this model, the resolution of the Gribov problem along 
the lines of Zwanziger's version of the stochastic approach jl] has been outlined. In the 
Zwanziger approach, the gauged configurations are attracted by the first Gribov regime. 
Here, it turned out that the approach can be modified in order to obtain attraction by the 
fundamental modular region. Moreover, a integrable drift force towards the fundamental 
region was developed. It was therefore possible to construct the gauge fixing action which 
generates the drift force in the Langevin simulation. The latter construction was essential 
to perform a heat bath simulation. In the case of the toy model, it was observed that the 
heat bath method is more efficient than the Langevin technique. 
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helpful discussions on gauge fixing and the stochastic approach. I am indebted to Oliver 
Schroder for valuable comments concerning, in particular, the construction of the conser- 
vative drift force. I thank R. Alkofer and M. Quandt for comments on the manuscript. 
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